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This study of the one dimensional Su-Schrieffer-Heeger model in a weak coupling perturbative 
regime points out the effective mass behavior as a function of the adiabatic parameter uj^/J, u)-n is 
the zone boundary phonon energy and J is the electron band hopping integral. Computation of low 
order diagrams shows that two phonons scattering processes become appreciable in the intermediate 
regime in which zone boundary phonons energetically compete with band electrons. Consistently, 
in the intermediate (and also moderately antiadiabatic) range the relevant mass renormalization 
signals the onset of a polaronic crossover whereas the electrons are essentially undressed in the 
fully adiabatic and antiadiabatic systems. The effective mass is roughly twice as much the bare 
band value in the intermediate regime while an abrupt increase (mainly related to the peculiar f D 
d ■ dispersion relations) is obtained at uj^ ~ \/2 J. 
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CIh' A sizeable electron-phonon interaction can induce a 
local deformation in the lattice accompanied by the for- 
mation of a quasiparticle with multiphononic character, 
the polaron H. As the spatial extension of the lattice 
deformation can vary, the concepts of large and small 
polaron have been introduced the transition be- 
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tween a large and a small polaron state is driven by the 
strength of the electron-phonon coupling and mon- 

itored through the behavior of ground state properties 
such as the polaron energy band and the effective mass 
Jl2| p"4[ . In particular, an abrupt increase of the polaron 
mass versus the e-ph coupling is associated with the oc- 
curence of the self-trapping event and the loss of the po- 
laron mobility properties. The aforementioned investiga- 
tions of polarons generally depart from the Hamiltonian 
of the Holstein model Jlq] in which the electron couples 
to dispersive optical phonons jl(| through a local (mo- 
mentum independent) interaction while the coupling to 
acoustical phonons, although possible in principle p7| , 
would lead to huge mass renormalizations [M. In fact, 
also the Holstein optical polaron masses are very heavy 
(at least larger than I0 3 times the bare band mass) in the 
^ , self-trapped state but masses of order 10 times the bare 
band mass are possible in the presence of high energy 
Q \ phonon spectra p9[ . 

As an alternative to the Holstein molecular crystal 
model with local interactions one may consider the Su- 
Schrieffer-Heeger (SSH) model Hamiltonian [fzoj . The 
tight binding SSH model was introduced to study poly- 
acetylene [ [2l| , an essentially one dimensional polymer 
with delocalized ir electrons which are responsible for 
the alternation of double and single bonds between two 
neighboring carbon atoms, a dimerized state p^ |. As a 
main characteristic of the SSH model Hamiltonian one 
notes that the electron-phonon interaction modifies the 
electron hopping matrix elements thus leading to a non- 
local (in momentum space) coupling with vertex function 
dependent on both the electron and the phonon wave 
vector. With the present paper we start an investigation 
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of the ground state and finite temperature properties of 
systems which can be described by the SSH model. Al- 
though some references to conducting polymers are made 
in the introductive notation we focus here on the general 
features of this model which may well apply to a class of 
quasi one dimensional systems with short range e-ph in- 
teractions p3|-p6|]. As a preliminar goal we evaluate the 
relevance of the mass renormalization (induced by this 
non local type of electron-phonon coupling) which might 
provide signatures of polaron formation in the system. 
In ID, the real space SSH Hamiltonian reads 

H = J2 - Wttf/r+l + fr+Jr) + 
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J nr+1 = --[j + a(u r -u r+1 j] (0.1) 

where J is the nearest neighbors hopping integral for 
an undistorted chain, a is the ir electron-phonon cou- 
pling, u r is the dimerization coordinate which specifies 
the displacement of the (CH) group on the r— lattice 
site along the molecular axis, p r is the momentum oper- 
ator conjugate to u r , M is the (CH) group mass, K is the 
effective spring constant, fj and f r create and destroy ir- 
electrons on the r— (CH) group. Incidentally we note 
that, by means of a Jordan- Wigner transformation, the 
Hamiltonian in eq.(l) maps onto that of a spin-Peierls 
chain [p7| , pq . After expanding the lattice displacement 
and its conjugate momentum in terms of the phonon cre- 
ation and annihilation operators b q and b q and Fourier 
transforming the electron operators, one gets the SSH 
Hamiltonian in momentum space: 

H = H + H int 

k q 

H int = E 9(k + q, k) (bl q + b q ) f\ +q f k 
k . q 

€k = —Jcos(k) 
J - 4^ sin 2 (I) 

9(k + q,k)= ^ 2MNuj {Mk + q)- sin(fc)) (0.2) 

where N is the total number of lattice sites. The 
phonon dispersion relation is defined in the range q € 
[0,7r]. Assuming a reduced Brillouin zone (\q\ < tt/2) 
the spectrum displays both an acoustic and an optical 
branch. Here we attack the Hamiltonian in eq.(2) by 
using the Matsubara Green's functions formalism and 
taking the e-ph term as the perturbation. Although a 
weak coupling perturbative approach cannot capture the 
full multiphononic nature of the polaronic quasiparticle, 
an envisaged, sizeable enhancement of the charge car- 
rier mass even to the lowest order of perturbative theory 
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would be an indicator of polaron formation. The full 
electron propagator is defined as: 

G(p,r) = - V(-l) n / dT 1 ...dT n (T T f p (T)H lnt (T 1 ) ■ 
n=o Jo X 

■H mt (r n )fl(0)) (0.3) 



where (3 is the inverse temperature, T T is the time 
ordering operator, < ... >o indicates that thermody- 
namic averages are taken with respect to the unperturbed 
Hamiltonian and only different connected diagrams con- 
tribute to any order n. 

I have calculated the self-energy terms due to one 
phonon (n = 2 in eq.(3)) and two phonons (n — 4 in 
eq.(3)) scattering processes. Their finite temperature ex- 
pressions are: 



n B {u q ) + n F (-e p ^ q ) 
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n B (oJ q ) + n F (e p - 
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g(p -q,p-q- Qx)gip -Qi,p-q-qi) 



G°„ ?1 (ie m - iui)G Q p _ q _ qi {ie rn - iu, 
Y, p 2h) (ie m ) = -^^2g 2 (p,P - q)g 2 (p - q,P - q - qi) ■ 
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D q {iLo n )D qi {iuoi) \G° p _ q {ie m - iu n ) 



4 2C) ( i£ m) = ^Ef 2 feP- ^9 2 (K k + q)> 
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^2 [ D q( iuJ n) Gp_ q (iE m - iU) n ) 

G° q+k ^e' m + iuj n )G° k (te' m ) 



E 



(0.4) 



There are three contributions due to different con- 
nected two-phonons diagrams p9[ . ns(w g ) is the Bose 
occupation factor and np(e p ) is the Fermi occupation 
factor. The frequencies s m and u) n are an odd and an 
even multiple of w/{3 respectively. G p (is m ) is the free 
electron propagator and D q (iuj n ) is the free phonon prop- 
agator. By analytical continuation ie rn — > e + iS one gets 
the retarded self-energy whose real part determines the 
renormalized electron mass rn e ff- 



m e ff _ 1 - <9-BeE p (e)/(9e|p = o ; E= _ 
m 1 + dReY. p (e)/de p \ p =o: e =- 



(0.5) 
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where, i?eE p (e) = ReT^\e) + ReY,f a \e) + 

ReY^p b \e) + ReY,^ c \e), has been obtained from eqs.(4) 
working out the straightforward but cumbersome dou- 
ble frequency summations and taking the zero tempera- 
ture limit. The model contains three free parameters: 
the hopping integral J, the zone boundary frequency 
oj v = 2y / K/M which coincides with the zone center op- 
tical frequency in the reduced zone scheme, the coupling 
constant a 2 /4K. In Fig.l, the mass ratio of eq.(5) is plot- 
ted versus the adiabaticity parameter w^/J assuming: (i) 
a narrow bare band J value, (ii) a weak e-ph coupling 

regime. n%h denotes the mass renormalization due to 

(2) 

the very one-phonon self-energy corrections while nT- e tf is 
enriched by the two-phonons scattering processes. Par- 
ticular care has to be taken in handling the principal val- 
ues which enter the real self-energy terms. I have used 
the representation 

RR(-)= lim^-T (0.6) 

and achieved numerical convergence by setting r; = 
10~ 4 and summing over 200 q points in each Brillouin 

zone. In the extreme adiabatic regime, ReT,fi\e) eval- 
uated at the band bottom is much larger than the two- 
phonons terms and the same trend holds for the par- 

(2) 

tial derivatives which enter the mass ratio. Then, Tn e ff 

does not differ essentially from m^ff and the renor- 
malization is very poor since scattering by low energy 
phonons does not enhance the electron mass over the bare 
band value. As an example, at lo^ — J/2, we get (in 
units meV) ReT^\-J) = 1.45, ReT^ a \~J) = 0.12, 
ReE^ h) (-J) = -0.17, Re^ c \-J) = -0.5- 10~ 3 . At 

2a) 

are comparable to ReT,^\— J) while the 2c diagram is 

(2) 

still negligible. Accordingly rn y e jfj starts to get larger 

than m^ff with an increase of ~ 15% in the interme- 
diate ~ J regime. Although, in this regime, the 
application of a low order perturbative theory may be 
questionable, the obtained mass enhancement is likely 
a signature that quasiparticlcs with polaronic character 
may form in the system once multiphonons scattering 
processes become a ppr eciable. Incidentally we note that 
a variational study |3(J] of the SSH model finds in the in- 
termediate regime favourable conditions for the existence 
of localized polaronic solutions. Fig. 2 displays the one- 
phonon effective mass versus to^/ J for two values of e-ph 
coupling. Setting J — 507716^ we take here an electron 
band narrower than in Fig.l. Again, the mass increase 
is relevant only in the intermediate Lu n /J- range while 
the electrons are essentially free in the extreme adiabatic 
and antiadiabatic regimes. The spike at ui^ ~ \/2J which 
dominates the mass behavior is mainly due to scattering 
by \q\ = 7r/2-phonons. This feature is related to the ID 
electron and phonon dispersion relations and it may be 



3 J/5, the ReT^ a \-J) and ReE^ h) (-J) terms 



4 



partly suppressed in higher dimensionality. In the ui n / J- 
window which is sensitive to renormalization effects the 
e-ph coupling parameter slightly affects the values of the 
one-phonon effective mass. For comparison also the one 
phonon effective mass of the Holstein-like model is re- 
ported on in Fig. 2. To point out the role of the e-ph 
coupling we have replaced the g 2 (p, p — q) function of the 
SSH model by a constant g 2 H without any change to the 
dispersion relations. Then the Holstein-like model here 
assumed differs from the "true" Holstein model jl5| of 
diatomic molecules in which purely optical phonons cou- 
ple to the electrons. We also note that simplified models 
with dispersionlcss Einstein phonon spectra would lead 
to wrong estimates of the ground state properties |l6|,0 ■ 
Setting g 2 H — a 2 /4K ■ uj^ with a 2 /4K = ImeV we ob- 
tain a mass behavior similar to that of the SSH model: 
again we find an abrupt mass increase at lo^ ~ y/2J (thus 
confirming its dependence on the ID dispersion relations 
rather than on the choice of the coupling function) but 
the present mass values are generally lower than those 
predicted by the Holstein model in strong coupling per- 
turbative theory . This suggests that the well known 
mass enhancement of the Holstein adiabatic and antia- 
diabatic polaron is due to the multiphononic effects fully 
captured (for instance) by the Lang-Firsov Q strong 
coupling method whereas the details of the short range 
coupling seem to have little influence on zone center prop- 
erties such as the effective mass. Conversely one might 
evaluate the SSH polaron mass in a strong coupling ap- 
proach to check whether and to which extent the mo- 
mentum dependent vertex function plays there a peculiar 
role. Finally, Fig. 3 emphasizes that the effective cou- 
pling a 2 /AK scarcely affect the mass renormalization in 
adiabatic and intermediate regimes while a slight mass 
dependence on a 2 /AK is observed in the antiadiabatic 
regime with very narrow electron band J = In 
the latter case however, in the upper portion of the x 
axis, the dimensionless coupling a 2 / (AK J) is larger than 
one and the perturbative method breaks down. In the 
intermediate case w T = J also the two-phonons mass is 
reported on. 

In conclusion, I have assumed a weak coupling per- 
turbative regime and applied low order diagrammatic 
techniques to the one dimensional Su-Schrieffer-Heeger 
Hamiltonian in order to compute the electron mass renor- 
malization induced by one- and two-phonons scattering. 
Tuning the parameter lo^/ J, the mass behavior has been 
analysed both in adiabatic, intermediate and antiadia- 
batic conditions. There is a sizable mass enhancement 
only in the intermediate and moderately antiadiabatic 
range whereas the electrons don't drag phonons whose 
energy is either much smaller or much larger than the 
electron energy. Hence, no mass enhancement is found 
in the extreme adiabatic and antiadiabatic regimes. Re- 
placing the momentum dependent vertex function by an 
Holstein-like coupling constant does not modify substan- 
tially the effective mass due to the one phonon self-energy 
diagram. We emphasize that this result holds in the 
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present weak coupling perturbative approach while differ- 
ent conclusions may be drawn in strong coupling theories. 
The two phonons self-energy corrections introduce an ap- 
preciable mass enhancement mostly in the intermediate 
u„ ~ J regime. This feature can be likely interpreted as 
an onset of polaronic crossover whose entity and precise 
location in parameter space requires however computa- 
tions of multiphonons effects. 



FIG. 1. Rcnormalized masses (in units of bare band elec- 
tron mass) versus the adiabaticity parameter, is due to 

(2) 

the one phonon self-energy correction. vn\^ includes also the 
two phonons self-energy terms. 
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FIG. 2. Rcnormalized masses (in units of bare band elec- 
tron mass) as due to the one phonon self-energy correc- 
tion versus the adiabaticity parameter. Two values of elec- 
tron-phonon coupling have been chosen. The mass behavior 
obtained in a Holstein-like model with constant coupling is 
reported on for comparison. 



FIG. 3. Mass renormalization (in units of bare band elec- 
tron mass) due to one phonon scattering versus the elec- 
tron-phonon coupling (in meV). Three values of the adia- 
baticity parameter have been chosen. In the intermediate 
case J — uj^ also the effect of the two phonons self-energy 
correction is displayed. 
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